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I report on recent progress and future prospects for lattice QCD calculations relevant for flavour physics and 
CP violation. I will focus on lattice studies that incorporate realistic vacuum polarization effects, i.e., with 
Uf = 2 + 1 sea quarks. 
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1. Introduction 

Theoretical calculations of nonperturbative QCD 
effects are essential to make full use of the re- 
sults generated by the experimental flavour program. 
The combination of theory and experiment can be 
used to extract the values of the elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix, look for 
new physics (NP) and constrain beyond the Stan- 
dard Model (BSM) theories. In order to be useful 
for achieving these goals, the theoretical calculations 
must have errors of the same order of the experimental 
errors. This requires uncertainties at the few percent 
level. 

Lattice QCD is a nonperturbative formulation of 
QCD based only on first principles. It also provides a 
quantitative calculation methodology, which has be- 
come a precise tool capable of providing some of the 
accurate determinations needed by phenomenology. 
The processes that can be analyzed with current lat- 
tice QCD methods in a precise way are those with 
stable (or almost stable) hadrons and no more than 
one hadron in the initial (final) state. Those processes 
include leptonic and semileptonic decays and neutral 
meson mixing. 

Accuracy in lattice calculations requires control 
over all the sources of systematic error. In particu- 
lar, it is essential to take into account vacuum po- 
larization effects in a realistic way, i.e., including up, 
down and strange sea quarks on the gauge configura- 
tions' generation. The up and down quarks are usu- 
ally taken to be degenerated, so those simulations are 
referred to as = 2 -I- 1. The vacuum polarization 
effects were almost always neglected in old lattice cal- 
culations due to limited computational power. This 
is known as the quenched approximation and intro- 
duces an uncontrolled and irreducible error, which can 
be as large as 10-30% [r|. Simulations with nj = 2 
sea quarks are still missing part of the vacuum po- 
larization effects and the associated systematic error 
is hard to estimate without repeating the calculation 
with Hf = 2 + 1 sea quarks. 

Another important source of systematic error is as- 
sociated with the fact that current simulations are un- 
able to simulate up and down quarks as light as the 
physical ones. The way of connecting lattice results to 
the physical world in a model independent way is by 



extrapolating those results using the guide of Chiral 
Perturbation Theory (ChPT). In order for this ex- 
trapolation to have controlled errors and be realistic, 
simulations must be performed for a range of light sea 
quark masses smaller than ms/2. 

Other systematic errors that must be included in 
any realistic lattice analysis are discretization, finite 
volume and renormalization effects. Discretization ef- 
fects can be estimated by power counting, but this 
estimate must be explicitly tested by performing the 
calculation at several values of the lattice spacing. 
Finite volume effects can be estimated by repeating 
the calculation at several volumes and/or using ChPT 
techniques. 

It is important to establish the validity of lattice 
methods by comparing its predictions against well 
known experimental quantities. Figure[l]shows the ra- 
tio between lattice results and experimental measure- 
ments for different decay constants, hadron masses 
and mass splittings. The agreement between lattice 
QCD and experiment is remarkable once vacuum po- 
larization effects are included in a realistic way. This 
gives us confidence in lattice techniques. 

In the next Sections I will discuss processes rele- 
vant for the experimental flavour physics program for 
which lattice QCD can provide accurate determina- 
tions of the nonperturbative inputs. I will thus re- 
strict my discussion to calculations with all sources of 
systematic error addressed. Among other things, that 
means that I will focus on simulations with Uf = 2 + 1 
sea quarks. 

1.1. Lattice fermion formulations 

The inclusion of quarks in the lattice QCD action, 
besides being expensive in computing time, has asso- 
ciated several difficulties. One of them is the so-called 
doubling problem, that consists of the fact that when 
one discretizes the QCD quark action in the most 
straightforward way, the naive discretization, there 
appear 15 additional unphysical tastes for any contin- 
uum flavour. Several methods exist to deal with the 
doubling problem. The most popular choice is (im- 
proved versions of) the Wilson action W\, that solves 
the doubling problem by adding a dimension five oper- 
ator with the cost of breaking chiral symmetry. Other 
fermion formulations keep an exact. Overlap or 
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Quenched with sea quarks 

Figure 1: Lattice QCD results divided by experimental re- 
sults in the quenched approximation (left panel) and with 
n/ — 2+1 sea quarks (right panel). Lattice values are 
obtained using improved staggered actions to describe sea 
quarks and u, d, s and c valence quarks, and NRQCD to 
simulate the 6 valence quarks. The plot is an updated 
version of the one in [31. 



almost exact, Domain Wall chiral symmetry pay- 
ing the price of significantly complicating the opera- 
tor structure and thus increasing the computing cost. 
The recent progress on the development of algorithms 
however, is starting to make feasible to perform un- 
quenched calculations with those formulations. 

Staggered fermions have good chiral properties and 
are computationally more efficient than any other 
light fermion formulation. The downside is that the 
doubling problem is not completely eliminated but re- 
duced to four tastes for any continuum flavour. In un- 
quenched simulations the extra tastes are eliminated 
by taking the root of the fermion determinant in the 
generation of configurations. This procedure is the fo- 
cus of intense scrutiny by the lattice community and 
there has been considerable progress in our under- 
standing of the issues involved during the past years 
. To date there is no proof of the correctness or in- 
correctness of the method, but all the tests performed 
on dynamical simulations taking the root of the deter- 
minant have given evidence of controlled effects which 
disappear in the continuum limit. 

The Asqtad (t} is the staggered formulation most 
widely used. It is improved with respect to the orig- 
inal staggered action so leading discretization effects 
are removed at tree-level. Remaining errors are there- 
fore 0{asa^) and 0{a'^). That is the fermion action 
that the MILC collaboration is using in the generation 
of the Hf — 2 + I configurations, which are employed 



in several of the lattice calculations referred to in the 
next sections. An even more improved staggered ac- 
tion is the recently developed HISQ (highly improved 
staggered quark) action Q . It reduces the remaining 
©(a^a^) errors coming from taste-changing effects in 
the Asqtad action by roughly a factor of three. 

For heavy quarks, charm and bottom, discretization 
errors coming in powers of the mass in lattice units, 
arriQ, are not negligible at typical lattice spacings. 

The HISQ action however incorporate corrections 
that remove the dominant amg effects so the lead- 
ing mass discretization effects are 0{as{amc)'^) and 
0{amc)'^. This action can thus be used in a very 
effective way to describe charm quarks if the lattice 
spacing is small enough^. 

To describe bottom quarks, however, an effective 
field theory framework is more adequate. Some im- 
plementations of the effective field methods are lattice 
heavy quark effective theory (HQET) (whose lead- 
ing term is the static approximation), non-relativistic 
QCD (NRQCD) and the Fermilab approach. The lat- 
tice NRQCD Lagrangian is obtained by discretizing 
the non-relativistic expansion of the continuum Dirac 
Lagrangian. The particular lattice NRQCD action 
used in the HPQCD calculations of B decay constants 
and -B*^ mixing parameters described in next sections 
is improved through 0{1/M^), 0{a'^) and leading rel- 
ativistic ©(l/M-^) The action parameters are fixed 
via heavy-heavy simulations, in particular the valence 
b quark mass is tuned to give the physical value of the 
T mass. 

The Fermilab action [HJ] is obtained following a dif- 
ferent approach. It starts with an improved relativis- 
tic Wilson action Q , which has the same heavy quark 
limit as QCD. With the Fermilab interpretation in 
terms of HQET \Td\ , this action can be used to describe 
heavy quarks without errors that grow as (amg)". 
One of the advantages of this approach is that it can 
be used for both charm and bottom quarks. The er- 
rors associated with the use of the Fermilab action 
mentioned in this paper, used by the FNAL/MILC 
collaboration, are O {asA/mq) and O (A/toq) . 



2. Leptonic decays 

The lattice determination of pseudoscalar decay 
constants, together with experimental measurements 
of pseudoscalar leptonic decay widths, can be used 
to extract the value of the CKM matrix elements in- 



^The MILC ensembles used by the HPQCD collaboration in 
its HISQ studies have small enough lattice spacings so the heavy 
discretization errors for the charm quark are under control. 
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volved in the process 

T{Pab ^lv) = (known factors) fp\Vab? ■ (1) 

On the other hand, for decay constants weU deter- 
mined experimentally, those for which the CKM ma- 
trix elements involved are known with a good pre- 
cision and experimental measurements are accurate, 
the comparison with lattice calculations can be used 
as a test of the theory used to make the theoretical 
prediction. 

The calculation of decay constants on the lattice is 
done using a simple matrix element, 



\qit,ibc\P{p)) = ifpPt, , 



(2) 



so they can be obtained with a very good accuracy. 
In the next subsections I will summarize the status of 
unquenched lattice calculations of charm and bottom 
decay constants. The current status of the lattice cal- 
culations of /k/ Itt, from which one can extract the 
value of the CKM matrix elementl T4s | , was reviewed 
by F. Mescia in this conference 



2.1. and Ds decay constants 

For the decays constants in the charm sector 
and there are lattice results available from two 
groups with nj = 2 + 1 sea quarks, the FNAL/MILC 
and the HPQCD 13] collaborations, to compare 
against experiment. Both use configurations gener- 
ated by the MILC collaboration for three different val- 
ues of the lattice spacing, a = 0.15 fm, a = 0.12 fm 
and a = 0.9 fm. The main difference between the two 
collaborations is the treatment of the valence quarks. 
While HPQCD uses the HISQ action for all the va- 
lence quarks, FNAL/MILC uses staggered Asqtad for 
the light quarks (up, down and strange) and the Fer- 
milab action for the charm quark. The HPQCD col- 
laboration has partially conserved currents, so they 
can extract the value of the decay constant without 
any renormalization. The FNAL/MILC collaboration 
needs to renormalize its currents, but they do it in 
a partially non-perturbative way that generates very 
small errors, around 1.5%. 

The values for /£>+ , fi)^ and the ratio of both quan- 
tities, together with the new CLEO-c experimental re- 
sults presented in this conference , are collected in 
Table H] Both lattice collaborations agree very well 
in the central values obtained for the decay constants. 
The errors of the HPQCD calculation are smaller than 
those for the FNAL/MILC due to the fact that the 
HISQ action is more improved^ than the Fermilab 



Group 


Reference 


/d+ 




fDjfD + 


experiment 


[14] 


{CLEO-c) 


(average) 


1.31(7) 


205.8(8.9) 


269.6(8.3) 


HPQCD 


[13J 


207(4) 


241(3) 


1.164(11) 


FNAL/MILC 


[12] 


215(14) 


254(14) 


1.188(26) 



Table I Comparison of /q+ and fo^ as obtained from ex- 
periment and from unquenched lattice QCD calculations. 



action and thus discretization errors are sensibly re- 
duced. The FNAL/MILC results agree with experi- 
ment for f£)+ and fo, within errors. The HPQCD 
result for /£>+ agrees also very well with the exper- 
imental number. However there is a discrepancy of 
over 3(7 between the HPQCD and experimental values 
for fu^ . This discrepancy has been recently suggested 
to be a. possible hint of beyond the Standard Model 
effects [15|. More work is needed to resolve this issue. 
From the experimental side, a reduction of the errors 
that will constrain the possible statistical fluctuations 
is expected. It is also desirable that certain issues like 
the assumption of three-generation CKM unitarity to 
set Vcs = Vud or the inclusion of radiative corrections 
from experimental data or Monte Carlo simulations 
are addressed. From the theory side, the error on 
the FNAL/MILC result for fo^ is currently larger by 
roughly a factor of four than the error on the HPQCD 
result. The FNAL/MILC collaboration plans to re- 
duce the uncertainties in their calculation in the near 
future by increasing statistics, using a smaller lattice 
spacing, and improving the determination of the in- 
puts needed. This will constitute an important check 
of the HPQCD numbers. 

Calculations of fD+ using Wilson fcrmions are also 
making progress, although they are still restricted to 
Nf = 2 simulations. A value of /d+ = 201 ± 22lg 
based on a single value of the lattice spacing was re- 
ported in p^ . The number is compatible with the 
Nf = 2-1-1 calculations but with larger errors. The 
authors in that reference determine the ratio /£>+ //^ 
and use the experimental value of /^r instead of di- 
rectly determine since the chiral logarithms are 
suppressed in the ratio and the associated system- 
atic error is thus reduced. Preliminary results using 
twisted mass fermions with Nf = 2 for fu+ and fn. 
were also reported by the ETM collaboration in [13] 
for two different values of the lattice spacing. 



^Improvement in this context refers to the addition of higher- 
dimensional operators to the action. 



2.2. B and Bs decay constants 

Lattice results for the decay constants in the B sec- 
tor are more needed than in the D sector since the 
corresponding CKM matrix elements to extract the 
information from experiment are worse known. The 
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value of the B decay constants are used in the SM 
predictions for processes very sensitive to beyond SM 
effects, such as Bg — > v^v~ . The purely leptonic de- 
cays themselves are also a sensitive probe of effects 
from charged Higgs bosons. 

The way of calculating these decay constants on the 
lattice is the same as for the charm decay constants. 
In fact, the FNAL/MILC collaboration has also cal- 
culated Jb , and the ratio of both with the same 
choice of actions, same ensembles and same procedure 
as for the D decay constants in ■ The errors in both 
the charm and bottom mesons are thus very similar. 
The results are listed in Table |TT1 In that table, the 
results from the other lattice rif = 2 + 1 calculation, 
by the HPQCD collaboration, are also included. In 
this case the HPQCD collaboration used the NRQCD 
action mentioned in section 11.11 to describe the b va- 
lence quark. Errors are then significantly larger than 
in their analysis of D'^ and Dg decay constants where 
they use the HISQ action for the c quark. A dominant 
source of uncertainty in their NRQCD calculation is 
the error associated with the one-loop renormalization 
applied. 



Group Reference fs Jb^ fss / fs 
FNAL/MILC fl2] 197(13) 240(12) 1.22(3) 
HPQCD [18] 216(22) 260(26) 1.20(3) 



Table II Values of /s and fos as obtained from the two 
lattice QCD calculations with n/ = 2 -I- 1. 

A suggested way of reducing the uncertainty in the 
calculation of Jb^ / Jb is by extracting it from the dou- 
ble ratio [fBs/fB]/[fK/fn] which can be calculated 
very accurately since it is very close to one in ChPT. 

3. Semileptonic decays 

Semileptonic decays can be used to extract CKM 
matrix elements like \Vcb\, |Kb|, \Vcd\, \Vcs\ and \Vus\- 
The theory input needed to get those parameters from 
experimentally measured semileptonic widths are the 
form factors in terms of which the hadronic matrix 
elements involved on those decays are parametrized. 
For example, for the decay D Klv, the differential 
decay rate is given by 

^ - (knownfactors)|K.|'/+(9') , (3) 

where f+{q^) is the vector form factor, which can be 
extracted from the matrix element of the vector cur- 
rent 

{K\V^\D) = U{q^){pD+PK^ + k{q^W . (4) 
with = (m|, — m'j()qf^ / q"^ . 



Lattice QCD can be used to calculate the value of 
those form factors as a function of the virtual W mo- 
mentum transfer, g^, or, equivalently, the recoil mo- 
mentum of the daughter meson. On the lattice, the 
smallest discretization errors correspond to the form 
factor at the largest momentum transfer, where the 
experimental data are less precise. In addition, the 
finite volume provides an infrared cutoff and there is 
a finite minimum value for the momenta that can be 
simulated. A way of circumventing this limitation is 
by using twisted boundary conditions that allow for 
arbitrary small values of the momenta [l^. Another 
set of important techniques that can be applied to 
semileptonic decay analysis are double ratio methods 
(20j . These methods can yield a reduction of both 
statistical and systematic uncertainties by a partial 
or total cancellation of those uncertainties between 
numerator and denominator in the ratios. 

In the case of |T4s |, experimental data for K nlv 
and lattice results for the corresponding form factors 
at zero momentum transfer yield a determination of 
this parameter with an error claimed to be 0.5 % -see 
F. Mescia's [ill and P. Gambino's [lH talks at this 
conference. 

The FNAL/MILC calculation of the shape of the 
form factor f'^ for the decay D — > Klv [24I in 2004 
constituted a prediction since the result appeared be- 
fore the experimental measurements by the FOCUS 
[23] and Belle ,24)] collaborations. The comparison of 
this lattice calculation with Belle experimental data 
is shown in Figure [5] as a function of q^ . The de- 
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Figure 2: Comparison of the shape of the vector form fac- 
tor f^(q'^) as measured experimentally and obtained on 
the lattice [25(|. 

termination of the vector form factor for D — > Klv 
and D — > nlv are being currently updated by the 
FNAL/MILC collaboration. The most straightfor- 
ward improvement is the extension of the calculation 
to smaller values of the lattice spacing. This will re- 
duce discretization errors which are the main source 
of uncertainty in the previous calculation. 
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The matrix element |Vcb| can be extracted from the 
decay B D*lv. The experimental results for this 
process at zero recoil have smaller errors than those for 
B Dlv. A new lattice calculation of the form fac- 
tor describing the decay, the axial vector form factor 
at zero recoil Ts^viXl, was presented in the lattice 
conference last year [2a]. This calculation eliminates 
the quenching errors from previous calculations since 
it includes n/ = 2 + 1 sea quarks. It also introduces 
a new double ratio method which gives the form fac- 
tor at zero recoil directly and with a reduction of the 
computational cost. The relation between the double 
ratio and the form factor, 



.,^|2_ {D*\nol^b\B){B\h,^^c\D*) 
'""'^ ^' {D*\nic\D*){B\hib\B) ' ^ ' 

is exact to all orders in the heavy-quark expansion 
in the continuum. Statistical errors in the numera- 
tor and denominator are highly correlated and largely 
cancel. And also, most of the renormalization can- 
cels, yielding a small uncertainty for the perturbative 
matching. 

The final result obtained for the form factor after 
chiral and continuum extrapolation is J^b_>£)*(1) = 
0.921 ± 0.013 ± 0.021 [23], where the first error is 
statistical and the second one includes all sources of 
systematic errors and is dominated by heavy-quark 
discretization errors. The CKM matrix element \Vcb\ 
extracted from this value of the form factor and the 
experimental averages in (28j is 



\Vcb\ = (39.2 ±0.6 ±1.0) X 10- 



(6) 



This value differs by 2(7 from the one extracted from 
inclusive decays. 

An alternative method for the extraction of \Vcb\ 
is the analysis of the decay B Dlv. A study of 
the form factors needed for such determination with 
quenched Wilson fermions was presented in [291 . An 
interesting aspect of that analysis is that the authors 
calculate the scalar form factor as well as the mo- 
mentum transfer dependence to avoid needing to ex- 
trapolate to zero recoil, where the experimental data 
suffer from phase space suppression compared to the 
B D*lv case. The scalar form factor, which only 
contributes for I = r can be used to constrain BSM 
physics [30j . 

The decay B irlv provides a way of extracting 
\Vub\ that is competitive with b ^ u inclusive decays. 
The only unquenched {Nf = 2 + 1) lattice determi- 
nation of the form factor needed for this extraction 
to date is the HPQCD calculation in [3l[, which uses 
the Asqtad staggered formulation for the light quarks 
and NRQCD for the b quarks. Together with exper- 
imental results, from this calculation one obtains 



\Vub\ = (3.55 ± 0.25 ± 0.50) X 10-3. The error is larger 
than it could be due to the fact that lattice and ex- 
perimental results for the corresponding decay rate 
have a poor overlap in q^. There are several meth- 
ods that can be adopted to try to address this issue. 
The HPQCD collaboration is reducing the recoil mo- 
mentum of the pion required for small by using a 
lattice frame in which the B meson is moving in the 
opposite direction to the pion. A modified version of 
NRQCD (moving-NRQCD) provides the description 
of the corresponding b quark with a large velocity 
[3^ . It is then possible to calculate the form factors at 
small q^ without needing large recoil momenta, hence 
keeping discretization and statistical errors under con- 
trol. The FNAL/MILC collaboration is implementing 
another method to overcome the poor overlap [33j . 
They are using the so called z— expansion to param- 
eterize the shape of the form factors. This param- 
eterization is model independent and based only on 
unitarity and analyticity [34| . When the experimental 
data are analyzed using the same model-independent 
parameterization, the results for the shape parameters 
from theory and experiment can be compared directly 
even if the respective q^ ranges have poor overlap. 
Hence, Lattice QCD must provide only the normal- 
ization of the form factor, which can be calculated 
at the q^ values where theoretical errors are smallest. 
The FNAL/MILC collaboration uses Fermilab action 
for the b quark and improved staggered formalism to 
describe the light quark in the numerical simulations 
for this work. The preliminary results obtained by 
this collaboration are shown in Figure [3l They look 
quite promising. 



3-5i I i I I I I I I I I 

— 3 param. fit constrained such that f^(0)=f^^(0) ~ x^/d.o.f. = 0.35 
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Figure 3: Preliminary results in [33| for the form factor 
describing the decay B — > irlu as a function of q^. 

The calculation of semileptonic and leptonic decays 
on the lattice can be used to construct ratios indepen- 



^ First error is experimental, the second one is the theoretical 



error from the lattice calculation of the vector form factor. 
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dent of CKM matrix elements, such 

T(D^^ii>) ' '^ith which one could test the consistency 
of lattice calculations against experiment or constrain 
BSM physics. 



4. Neutral meson mixing 

The theoretical determination of the parameters 
that describe the mixing in the neutral Kaon and B 
systems are needed in order to test for the consistency 
of the SM description of CP violation when compared 
against experiment. 

In the SM, the neutral K and B meson mixing is 
due to box diagrams with exchange of two M^-bosons. 
These box diagrams can be rewritten in terms of an 
effective Hamiltonian with four-fermion operators de- 
scribing processes with AF — 2 (F=S or B). The 
matrix elements of the operators between the neutral 
meson and antimeson encode the non-perturbative in- 
formation on the mixing and can be calculated using 
lattice QCD techniques. 

4.1. Indirect CP violation in neutral kaon 
decays: Bk- 

The non-perturbative input to study CP violating 
effects in K'^ — mixing is parametrized by i?^, 
defined as 



(j^°|QA5^2(M)|j^°) 

f(^°|s7M75d|0)(0|s-7^75d|^°> 



(7) 



The theoretical calculation of this parameter, to- 
gether with the experimental measurement of sk = 



A(gs^(7r7r)j_o) ' S^^^^ ^ hyperbole in the p - 77 plane, 
where p and 77 are the usual unitarity triangle pa- 
rameters. That corresponds to the light blue band in 
Figure [S] 

The most recent unquenched lattice determinations 
of i3if(/x) are shown in Figured The current most 
accurate value is the one by the RBC/UKQCD collab- 
oration using domain wall fermions with n/ = 2 + 1 
sea quarks [37] 

5^(2 GeV) = 0.524 ± 0.010 ± 0.013 ± 0.025 , (8) 

where the first error is statistical, the second one is 
associated with the non-perturbative renormalization 
and the third one corresponds to the other systematic 
errors. The main source of uncertainty in this cal- 
culation is discretization errors. The result in ([5]) is 
obtained from simulations with a single value of the 
lattice spacing, = 1.729(28)GeV. The discretiza- 
tion error corresponding to the use of this single value 
of a is estimated in [s^l by using the scaling behaviour 
of a previous quenched calculation with the same light 



' 1 ' 1 
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Hporn/T iKOrn 7mf, l«fa.>f. l 




Nj = 2+1, lenormalization P 




RBC/UQCD 2007 (DW) 





= 2+1, renonnalization NP 





Figure 4: Recent unquenched lattice values of Bk- The 
parameter is given in the MS — NDR scheme and 
at an scale equal to 2GeV. The results are taken 
from [111 (JLQCD), Q (HPQCD/UKQCD) and ^ 
(RBC/UKQCD). The matching is done nonperturbatively 
in [3^ and [s^l, and perturbatively at one-loop in [s^. The 
perturbative matching is the origin of the rather large er- 
ror. 



quark action. The authors found that error to be of 
around a 4%, which dominates the systematic uncer- 
tainty. The RBC/UKQCD collaboration is planning 
to explicitly check discretization effects by performing 
the unquenched calculation with the same quark and 
giuon actions but with a smaller lattice spacing. 

The JLQCD result in Figure H with iV/ = 2 sea 
quarks is affected by the fact that part of the vacuum 
polarization effects are still missing and are difficult to 
estimate a priori. This collaboration is planning on 
extending its methodology to a Nj = 2-1-1 calculation 
of Bk- The needed configurations are being currently 
generated. The ensemble will include also configura- 
tions for two different volumes. This will allow them 
to explicitly study and reduce the finite volume ef- 
fects that are the main source of uncertainty in their 
Nf = 2 calculation. 

There are two other unquenched determinations of 
Bk in progress which use mixed actions*. A Nf — 2 
calculation outlined in 38] uses overlap fermions for 
the valence quarks and twisted mass fermions for the 
sea quarks. A Nf =2-1-1 study whose preliminary 
results can be found in [s^ uses domain wall valence 
fermions and improved staggered sea fermions. The 
errors due to the matching to the continuum affect- 
ing previous staggered calculations can be highly re- 



*In a mixed action calculation the valence and sea quarks 
are described with different fermion formulations. 
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duced in this analysis by performing the renormal- 
ization non-perturbatively. The use of domain wall 
valence quarks makes the chiral extrapolation more 
continuum-like than in the purely staggered case. An- 
other advantage of this calculation is that there are 
staggered configurations generated by the MILC col- 
laboration for a large range of lattice spacings, vol- 
umes and small sea quark masses. This allows for 
good control over the systematic error from chiral and 
continuum extrapolations. The expected final error 
from this calculation is around 5%. 

In the near future, there will be thus several lattice 
calculations of using different discretizations and 
with errors at the 5% level. 



4.2. 5° mixing: AM^,^, Ar^,^ and ^ 

The mixing in the — B^ system is an interest- 
ing place to look for NP effects. The BSM effects can 
appear as new tree level contributions, or through the 
presence of new particles in the box diagrams. In fact, 
it has been recently claimed that there is a disagree- 
ment between direct experimental measurement of the 
phase of B^ mixing amplitude and the SM prediction 
[40| . Possible NP effects have also been reported to 
show up in the comparison between direct experimen- 
tal measurements of sin(2/3) and SM predictions us- 
ing i?° mixing parameters [4l|. Studies of neutral B 
meson mixing parameters can also impose important 
constraints on different NP scenarios 42]. 

The quantities that describe the mixing in the i?^ 
system are the mass differences, AMg^^, and decay 
width differences, AFs^d, between the heavy and light 
and mass eigenstates. The non-perturbative 
physics of those processes is contained in hadronic 
matrix elements of the four-fermion operators in the 
effective Hamiltonian with AB — 2. Those matrix el- 
ements are parametrized by products of B decay con- 
stants and bag parameters. For example, for the mass 
difference 



AM, 



's(d)\theor. OC | Vj* (^) Vtfc ^/j 



Bf 



(9) 



with (50|Qf)|i?, 



M 



s{d} 



(d) 



and 



O 

Many of the uncertainties that affect the theoreti- 
cal calculation of the decay constants and bag param- 
eters cancel totally or partially if one takes the ratio 

= BbJJb^Bbj^- Hence, this ratio and therefore 
the combination of CKM matrix elements related to it, 
I I , can be determined with a significantly smaller 
error than the individual matrix elements. The ra- 
tio ^ is also an important ingredient in the unitarity 
triangle analyses. 

The first lattice calculation of the B'^ mixing param- 
eters with Hf = 2-1-1 sea quarks, which only studied 



the sector, was performed by the HPQCD collab- 
oration in ^S]. The authors obtained 

AMs = 20.3(3. 0)(0.8)ps"i and AF, = 0.10(3)ps-i . 

(10) 

which is compatible with experiment. 

The FNAL/MILC and HPQCD collaborations are 
currently working on a more complete study of B'^ — 
B^ mixing, including B^ and B^ parameters. The 
main goal of both projects is obtaining the ratio ^ fully 
incorporating vacuum polarization effects. The choice 
of actions and the setup is the same as for their fs, 
Jb^ calculations described in Section [2] -more details 
can be found in [3]. 

Results from the two collaborations are still prelim- 
inary. Figures [5] and [S] show some examples of the 
values of /bV^^bBb obtained as function of the light 
sea quark masses. The dependency on the light sea 



0.9 
0.85 
> 0.8 

o 

0.75 

f 

ta" 0.65 
5- 0.6 
055 
0.5, 



Preliminary 



^ a— 0. 1 2 fm (improving .stati.stics) 
■ a=0.12fm 
▲ a=0.09fm 



0.4 

/m 



Figure 5: Values of /s, \/Mb^BbI in GeV'^''^ as a function 
of the light sea quark mass normalized to the physical 
strange quark mass from the HPQCD collaboration. The 
data include statistical, perturbative and scale errors. The 
bottom valence quark is fixed to its physical value and the 
strange valence quark is very close to its physical value. 
The strange sea quark mass is also very close to its physical 
value. 

quark mass is in both studies very mild, so only the 
chiral extrapolation in the d quark mass for B^ param- 
eters is expected to be a significant source of error. In 
Figure [5l it can also be appreciated that the results 
for the two different lattice spacings are very similar, 
which indicates small discretization errors. 

A comparison of the preliminary results from the 
two collaborations for the ratio ^ is shown in Figure [T] 
Only the full QCD points^ are included. The results 
of the two collaborations agree within statistical er- 
rors. This is very encouraging since both analyses use 



^Where valence and sea quark masses are the same. 
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■ 
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"=0.007/0.050 


■ 
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Preliminary 



Figure 6: Bare values of /s, Mb^Bb^ in lattice units as 
a function of the light sea quark mass also in lattice units 
from the FNAL/MILC collaboration. Results are shown 
for both and including only statistical errors. The 
results correspond to one of the three lattice spacings at 
which the FNAL/MILC's study is performed. The bot- 
tom valence quark is fixed to its physical value and the 
strange valence quark is very close to its physical value. 
The strange sea quark mass is also very close to its physi- 
cal value. 



The effects of heavy nevif particles in the box di- 
agrams that describe the mixing can be seen in 
the form of effective operators built with SM degrees 
of freedom. The NP could modify the Wilson coeffi- 
cients of the four-fermion operators that already con- 
tribute to -B" mixing in the SM and gives rise to new 
four-fermion operators in the A_B — 2 effective Hamil- 
tonian -see [3, Ell ^ ^i^^ ^^'^ possible operators 
in the SUSY basis. The calculation of those Wilson 
coefficients for a particular BSM theory, together with 
the lattice calculation of the matrix elements of all the 
possible four-fermion operators in the SM and beyond 
and experimental measurements of mixing param- 
eters, can constraint BSM parameters and help to un- 
derstand new physics. To date, there does not exist an 
unquenched determination of the complete set of ma- 
trix elements of four-fermion operators in that general 
AS = 2 effective Hamiltonian. However, the two col- 
laborations currently working on mixing in the SM 
are planning to extend their analysis to BSM opera- 
tors in the near future. Actually, the HPQCD collab- 
oration has already calculated the one-loop matching 
coefficients needed for such an analysis |47| . 



5. Conclusions 



Preliminary 



■ a=O.I2fm, HPQCD 

▲ a=0.09 fm, HPQCD 

a=0. 1 2 fm, FNAL/MILC 

y a=0.09 fm, FNAL/MILC 



Figure 7: Product of the ratios ^ and Mb, /Mb^ as a func- 
tion of the down quark mass normalized to the strange 
quark mass. Resuhs for both FNAL/MILC and HPQCD 
collaborations including only statistical errors are shown 
for two different values of the lattice spacing a. 



completely different descriptions for the heavy quarks. 

The final step in those calculations is extrapolating 
the results to the physical values of the light quark 
masses and, in the case of the FNAL/MILC collabora- 
tion, performing simultaneously the extrapolation to 
the continuum. Systematic errors on those extrapola- 
tions are currently being studied. These two analyses 
are expected to have final results very soon with total 
errors ranging 5 — 7% for Jb^ \/^b, and 2 — 3% for 



Hints of discrepancies between SM predictions and 
experimental measurements have started to show up 
in some CP violating observables [1^. As claimed 
in (48j . the precise determination of parameters like 
Bk, Ik and ^, and CKM matrix elements like \ Vcb\ is 
crucial in order to fully exploit the potential of CP vi- 
olating observables on constraining NP. Lattice QCD 
has a fundamental role in that program. There has 
recently been important progress in order to achieve 
realistic lattice calculations, with n/ = 2-1-1 sea quarks 
and a serious study of systematic errors. New results 
relevant for phenomenology have appeared in the last 
year in the Kaon and D meson sectors with errors at 
the few percent level. In the near future, results for 
i?" mixing parameters and B decay constants will be 
also available with errors at the few percent level. 

Figure [5] summarizes the impact of the lattice cal- 
culations with Uf 2 + 1 on the unitarity triangle 
analysis. In generate the plot, the value of Bk is the 
one by the RBC/UKQCD collaboration [H, is 
taken from leptonic decays with Jk/ fir given by the 
HPQCD collaboration [13|, \Vcb\ is from semileptonic 
B ^ D*lv with the form factor by the FNAL/MILC 
collaboration [27*1, \Vub\ is from Flynn and Nieves [49| 
using, among other information, the form factor for 
B ttIv by the HPQCD ^ and FNAL/MILC 
[H, [lOl collaborations. A final Uf = 2-1-1 result for ^ 
is not yet available. In order to illustrate its effect on 
the p — rj plane in Figure 8, we assumed a value for ^ 
with a 3% error. As explained in Section [321 this is 
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Figure 8: Constraints on the p — f) plane imposed by recent 
unquenched lattice QCD calculations. The black point 
and error bars correspond to the values of p and f] from 
the PDG 2006 review of particle physics [Sll • See the text 
for more explanations. The plot is an updated version of 
the one in 



the expected error for ^ from both the HPQCD and 
the FNAL/MILC calculations. 

Several lattice collaborations are currently produc- 
ing accurate nj = 2 + 1 results as discussed in this 
paper, and other collaborations are starting to gener- 
ate n/ = 2-1-1 ensembles using different sea quark ac- 
tions. This will allow an important consistency check 
of lattice methods. 
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